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Abstract 

The magnetic structure of LaMn03 is investigated on three-dimensional clusters of 
MnOg octahedra by using a combination of relaxation and Monte Carlo techniques. It is 
found that the cooperative Jahn- Teller phonons lead to the stabilization of A-type anti- 
ferromagnetic and C-type orbital structures in the physically relevant region of parameter 
space for LaMn03 with small corrections due to tilting effects. The results suggest that 
strong Coulomb interactions are not necessary for a qualitative description of undoped 
manganites. In fact, it is shown that the present result is not essentially changed even if 
the Coulomb interaction is explicitly included. 

1 Introduction 

The study of manganese oxides is receiving considerable attention in recent years both in its 
theoretical and experimental aspects M. From the technological viewpoint, these materials 
could be used in the preparation of high-sensitive magnetic-field sensors due to their colossal 
magnetoresistance (CMR) phenomena. In addition, researchers in the condensed matter 
field have been interested in the rich phase diagram of these materials originating from the 
competition and interplay among charge, spin, and orbital degrees of freedoms. Obtaining a 
unified picture for this rich phase diagram is a challenging open problem. 

A prototype for the theoretical investigation of manganese oxides is the double-exchange 
(DE) framework, describing the hopping motion of e g -electrons ferromagnetically coupled to 
localized t2 g -spins. This idea has conceptually explained the appearance of ferromagnetism 
when holes are doped ||. In addition, within the one orbital model, the existence of phase 
separation has been recently unveiled with the use of modern numerical techniques ||, leading 
to a potential explanation of the CMR effect Q . 

However, in order to understand the fine details of the phase diagram of manganites, 
the one-orbital model is not sufficient since the highly nontrivial A-type spin antiferro (AF) 
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and C-type orbital structures observed experimentally in the undoped material LaMn03^] 
cannot be properly addressed in such a simple context. Certainly two-orbital models are 
needed to consider the nontrivial state of undoped manganites. In this framework the two- 
band model without phonons has been studied before, and the importance of the strong 
Coulomb repulsion has been remarked for the appearance of the A-AF state [§, 0, |, g. This 
is based upon the belief that the competition between kinetic and strong correlation effect 
determines the optimal orbital for e g -electrons and the lattice will be simply distorted to 
reproduce such optimal orbitals. However, Coulombic approaches have presented conflicting 
results regarding the orbital order that coexists with the A-type spin state, with several 
approaches predicting G-type orbital order, which is not observed in practice. 

While it is certainly correct that the orbital degrees of freedom play an essential role for 
the stabilization of A-AF, it should be noticed that our understanding is still incomplete. 
In particular, it is important to remark that the orbital structure is tightly related to the 
Jahn- Teller (JT) distortion of MnOe octahedron. If each JT distortion would occur indepen- 
dently, optimal orbitals can be determined by minimizing the kinetic and interaction energy 
of e g -electrons. However, oxygens are shared between adjacent MnOg octahedra, indicating 
that the JT distortions occurs cooperatively. Especially in the undoped situation, all MnOg 
octahedra exhibit JT distortion, indicating that the cooperative effect is very important, as 
discussed by Kanamori [10|. Thus, in order to understand the magnetic and orbital structures 
in LaMn03, it is indispensable to optimize simultaneously the electron and lattice systems. 
However, not much effort has been devoted to the microscopic treatment of the cooperative 
effect 1 11], although the JT effect in the manganese oxide has been studied by several groups 
12, |l^, 14, 15, 16]. Then, in the present work, a careful investigation of this problem is per- 



formed with some numerical techniques, focusing on n = 1, where n is the electron number 
per site. 

In this paper, the optimal oxygen positions are determined by a relaxation technique 
to obtain the lattice distortions corresponding to several t2 g -spin magnetic structures. In 
addition, Monte Carlo (MC) simulations were also performed to investigate the spin and 
orbital structure without a priori assumptions for their order. It is found that A-AF, as well 
as the C-type orbital structure, occurs in realistic parameter regions for LaMn03, i.e., large 
Hund coupling between the e g -electron and t2 g -spin, small AF interaction between t2 g -spins, 
and strong electron-lattice coupling. It should be emphasized that our results are obtained 
without the Coulomb interaction. It is shown in a simple case that the optimized results are 
essentially unchanged even if the Coulomb interaction is included explicitly in the model. 

The organization of this paper is as follows. Section 2 is devoted to the formulation to 
include the cooperative effect in the two-orbital model tightly coupled to the JT distortion, 
and some technical points are briefly discussed. In Sec. 3, the results on the magnetic and 
orbital structures are provided and it is shown that the region of A-AF in the magnetic phase 
diagram is reasonable for LaMn03, since the couplings needed agree with experiments. In 
Sec. 4, a prescription to obtain the C-type orbital order with the alternation of 3x 2 — r 2 and 
3y 2 — r 2 orbitals is provided. Finally, in Sec. 5, the effect of the Coulomb interaction on the 
lattice distortion is discussed in the ferromagnetic state. Throughout this paper, units such 
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that h = k-Q = 1 are used. 



2 Formulation 
2.1 Hamiltonian 

Let us consider the motion of e g -electrons tightly coupled to the localized t2 g -spins and the 
local distortions of the MnC>6 octahedra. This situation is well described by 

H = H 2 oTb + -ffAFM + #el-ph + -^el-el- (1) 

Here the first term indicates the two-orbital Hamiltonian, given by, 

#2orb = - ^2 i 77' c i 7 <r C i+a7' CT ~ <^H ^2 Si ' C \-ya a ™' C n°' > ( 2 ) 
ia77'cr i-yaa' 

where c- ma (cu, a ) is the annihilation operator for an e g -electron with spin a in the d x 2_ y 2 
(d 5z 2_ r 2) orbital at site i. The vector connecting nearest-neighbor sites is a, t* / is the 
hopping amplitude between 7- and 7'-orbitals connecting nearest-neighbor sites along the a- 
direction via the oxygen 2p-orbital, Jh is the Hund coupling, Si the localized classical t2 g -spin 
normalized to |Sj| = 1, and a = (a±, (T3) are the Pauli matrices. 

The second term is needed to account for the AFM character of the manganese oxide, 
given by 

-Hafm = J' ^ Sj • Sj, (3) 
<iJ) 

where J' is the AF-coupling between nearest-neighbor t2 g -spins. 

In the third term, the coupling of e g -electrons to the distortion of MnC>6 octahedron is 
considered as 

#el-ph = 9 X! ^(Qli^O + ^2i^l + fti^JiyCif a 
i<T 77 ' 

+ (l/2)^[fc br Q? i + A;jT(Q| i + Qi i )], (4) 

i 

where g is the electron-phonon coupling constant, Qn denotes the distortion for the breathing 
mode of the MnOe octahedron, Q 2 i and Qs\ are, respectively, JT distortions for the (x 2 — y 2 )- 
and (3z 2 — r 2 )-type modes, and o"o is the 2x2 unit matrix. Spring constants for breathing- 
and JT-modes are denoted by febr and /cjt, respectively. 

The final term indicates the Coulomb interactions between e g -electrons. As mentioned in 
Sec. 1, since only the JT phonons are investigated, the Coulomb interactions are neglected. 
This point will be discussed later in the text (Sec. 5) where the explicit form of i? e i-ei is 
provided and briefly analyzed. 
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Figure 1: (a) 2 x 2 x 2 lattice composed of eight Mn06 octahedra. Note that oxygens are 
shared by adjacent octahedra. (b) Three kinds of distortions of octahedron considered in this 
paper. The arrows indicate the direction of displacement of oxygens. 

2.2 Lattice distortion 

As shown in Fig. 1(a), oxygens are shared between adjacent octahedra, indicating that the 
local lattice distortions cannot be treated independently and a cooperative analysis is needed 
for this problem. For this purpose, the normal coordinates for distortions of the MnOg 



octahedron, shown in Fig. 1(b), are written as [11] 



for the breathing mode, 
and 



Qu = (l/v / 3)(A xi + A yi + A zi ), (5) 
Q 2i = (l/V2)(A xi - A yi ), (6) 



Q 3i = (l/^6)(2A zi -A xi -A yi ), (7) 
for the JT modes, where A a j is given by 

A ai = A a + J ai . (8) 

The first term indicates the deviation from the cubic lattice, given by A a = L a — L, where L a 
is the length between Mn-ions along the a-axis and L = (L x + L y + L z )/3. The second term 
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is the contribution from the shift of oxygen position, expressed by <5 a i = uf — uf , where 
uf is the deviation of oxygen from the equilibrium position along the Mn-Mn bond in the 
a-direction. By this consideration, the cooperative JT distortion as well as the macroscopic 
lattice deformation is reasonably taken into account. Note that the buckling and rotational 
modes of MnOg octahedron are not explicitly included in this work. In general, L a can be 
different for each direction, depending on the bulk properties of the lattice. Since the present 
work focuses on the microscopic mechanism for A-AF formation in LaMn03, the undistorted 
lattice with L x = L y = L z is treated first, and then corrections will be added. 

2.3 Hopping amplitudes and energy scale 



In the cubic undistorted lattice, the hopping amplitudes are given by [17| 



for the x-direction, 
for the y-direction, and 



& = - V3£ b = = 3£ b = t, (9) 

& = ^ = y^L = 3t bb = t, (io) 

*bb = 4t/3, = *ab = *ba = 0> (H) 



for the z-direction. Throughout this paper, the energy unit is t. 

Corresponding to this choice of the energy unit, the length in the lattice distortion is 
scaled by y/t/kjT- As a result of this scaling, a non-dimensional electron-phonon coupling 
constant A is defined as 

A = g/Vk^i. (12) 

It is noted that this coupling constant can be related to the static JT energy, which is 
conventionally defined by Ejt = g 2 /(2kj^), as 

£ JT = t\ 2 /2. (13) 

Note also that the present length scale is rewritten as 

\ft/k^ = WA, (14) 

where £jt = g/kjT is the characteristic length for the JT distortion. From the experimental 
result, £jt is estimated as 0.3A ||, which is a typical length in this context. 

As for the spring constant for the breathing mode, it is expressed as fcbr = /%jt and 
the ratio (5 is treated as a parameter. If it is plausibly assumed that the reduced masses 
for those modes are equal, this ratio is given by f3 = (lo^/lojt) 2 , where Wbr and wjt are the 
vibration energies for manganite breathing- and JT-modes, respectively. From experimental 
results and band-calculation data, Wbr and wjt are, respectively, estimated as ~ 700cm _1 and 
500-600cm _1 [p^j . Then, throughout this work, (3 is taken as 2, although the results presented 
here are basically unchanged as long as (3 is larger than unity. 
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Figure 2: Magnetic structures in 2 x 2 x 2 clusters. 

In this work, the change in t due to the displacement of oxygen position is not taken 
into account, but such an effect is shown to be very small as follows. Due to the pseudo- 
potential theory fl"9|| , the exact hopping amplitude, for example, along the x-direction between 
a-orbitals in i and i + x sites, is expressed as 



r 



t 



(1-e 2 ) 7 / 2 ' 



(15) 



with e = \2uf\/L x . It should be remarked that the change in due to the oxygen shift is 
of the order of e 2 , not of the order of e. Since e is estimated to be at most a few percent, 
the change is considered to be negligible. Thus, in the present work, such a change is not 
included to avoid unnecessary complication in the calculation. However, when the distorted 
lattice is considered, namely, when the deviation from L x = L y = L z is taken into account, 
the change in the hopping matrix due to this distortion should be included, because the effect 
is of the order of A a /L in this case, not of the order of (A a /L) 2 This point will be discussed 
again in Sec. 4. 



2.4 Techniques 

To study Hamiltonian Eq. (||) without H e \_ e \, two numerical techniques have been applied. 
One is the relaxation technique, in which the optimal positions of the oxygens are determined 
by minimizing the total energy. In this calculation, only the stretching mode for the octahe- 
dron, is taken into account. Moreover, the relaxation has been performed for fixed structures 
of the t 2g -spins such as ferro (F), A- type AF (A-AF), C-type AF (C-AF), and G-type AF 
(G-AF), shown in Fig. 2. The advantage of this method is that the optimal orbital struc- 
ture can be rapidly obtained on small clusters. However, the assumptions involved in the 
relaxation procedure should be checked with an independent method. 

Such a check is performed with the unbiased MC simulations used before for one- and 



two-dimensional clusters using non-cooperative JT-phonons [15]. The dominant magnetic and 
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Figure 3: Total energy as a function of J' on a 2 x 2 x 2 lattice with Jh = 8 and A = 1.5. The 
solid lines and circles indicate the relaxation and MC results, respectively. MC simulations 
have been performed at temperature 1/200. 

orbital structures are deduced from correlation functions. In the MC method, the clusters 
currently reachable are 2x2x2, 4x4x2, and 4x4x4. In spite of this size limitation, 
arising from the large number of degrees of freedom in the problem, the available clusters are 
sufficient for our mostly qualitative purposes. In addition, the remarkable agreement between 
MC and relaxation methods lead us to believe that our results are representative of the bulk 
limit. 

3 Results 

3.1 Magnetic structure 

In Fig. 3, the mean-energy is presented as a function of J' for Jh = 8 and A = 1.5, on a 2 x 2 x 2 
cluster with open boundary conditions. The solid lines and symbols indicate the results 
obtained with the relaxation technique and MC simulations, respectively. The agreement is 
excellent, showing that the relaxation method is accurate. The small deviations between the 
results of the two techniques are caused by temperature effects. As intuitively expected, with 
increasing J' the optimal magnetic structure changes from ferro- to antiferromagnetic, and 
this occurs in the order F^A-AF^C-AF^G-AF. 

To check size effects, the t2 g -spin correlation function S(q) was calculated also in 4 x 4 x 2 
and 4x4x4 clusters, where 




(16) 



Here N is the number of sites and (• • •) indicates the thermal average value. As shown in 
Fig. 4, with increasing J' the dominant correlation changes in the order of q = (0,0,0), 
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Figure 4: Spin correlation function 5(q) obtained by MC simulations as a function of 
J', at Jh = 8 and A = 1.5. Solid and open symbols denote the results in 4 x 4 x 2 and 
4x4x4 clusters, respectively. Circles, squares, diamonds, and triangles indicates 5(q) for 
q= (0,0,0), (tt, 0, 0), (tt, it, 0), and (7r,7r,7r), respectively. 

(vr,0, 0), (jt, 7r, 0), and (tt, it, it). The values of J' at which the spin structures changes agree 
well with those in Fig. 3. 

3.2 Orbital structure 

The shapes of the occupied orbital arrangement with the lowest energy for F, A-AF, C- 
AF, and G-AF magnetic structures are shown in Fig. 5. For the F-case, the G-type orbital 
structure is naively expected, because it is believed that the ferromagnetic spin structure is 
favored by the AF orbital configuration. However, a more complicated orbital structure is 
stabilized in the actual calculation, indicating the importance of the cooperative treatment 
for JT-phonons. For the A-AF state, only the C-type structure is depicted in Fig. 5, but 
the G-type structure, obtained by a 7r/2-rotation of the upper x-y plane of the C-type state, 
was found to have exactly the same energy. Small corrections will remove this degeneracy in 
favor of the C-type as described in the next section. For C- and G-AF, the obtained orbital 
structures are G- and C-types, respectively. 

Note that there exists an additional triplet degeneracy in the A-type state due to the cubic 
symmetry for each magnetic structure: If axes are changed cyclically (x^y,y^z,z^x), 
the optimized orbital structure is also transformed by this cyclic change, but the energy is 
invariant. Then, the magnetic and orbital structure in LaMnOs occurs through a spontaneous 
symmetry breaking process. 

Although the same change of the magnetic structure due to J' was already reported in 
the electronic model with purely Coulomb interactions Q, the orbital structures in those 
previous calculations were G-, G-, A-, and A-type for the F-, A-AF, C-AF, and G-AF spin 
states, respectively. Note that for the A-AF state, of relevance for the undoped manganites, 
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Figure 5: Optimized orbital structure for each magnetic structure. 

the G-type order was obtained ||, although in another treatment for the Coulomb interaction, 
the C- and G-type structures were found to be degenerate |7]], as in our calculation. Thus, 
the stabilization in experiments of the C-type orbital structure is still puzzling both in the 
JT and Coulomb mechanisms. This point will be discussed later in the text. 

3.3 Magnetic phase diagram 

In Figs. 6(a) and (b), the phase diagrams on the (J',A)-plane are shown for Jh = 4 and 
8, respectively. The curves are drawn by the relaxation method. As expected, the F-region 
becomes wider with increasing Jh- When A is increased at fixed Jh, the magnetic structure 
changes from F^A-AF— >C-AF^G-AF. This tendency is qualitatively understood if the two- 
site problem is considered in the limit Jh > 1 and i?jT 3> 1. The energy-gain due to the 
second-order hopping process of e g -electrons is roughly 5E\p ~ 1/Jh and SEp ~ 1/-Ejt 
for AF- and F-spin pairs, respectively. Increasing Ejt, 5Ep decreases, indicating the relative 
stabilization of the AF-phase. In our phase diagram, the A-AF phase appears for A > 1.1 and 
J' < 0.15. This region does not depend much on Jh, as long as Jh 3> 1. Although A seems 
to be large, it is realistic from an experimental viewpoint: Ejt is 0.25eV from photoemission 
experiments [2C] and t is estimated as 0.2 ~ 0.5eV [21], leading to 1 < A < 1.6. As for J', it 
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Figure 6: Magnetic phase diagram on the (J', A) plane for (a) Jh = 4 and (b) 8 (relaxation 
method). Below the thin solid lines in the F and A-AF regions, the JT-distortion disappears, 
suggesting that the system becomes metallic. 

is estimated as 0.02 < J' < 0.1 (||, [22| . Thus, the location in parameter-space of the A-AF 
state found here is reasonable when compared with experimental results for LaMnOs. 

4 Orbital order in the A-AF phase 

Let us now focus on the orbital structure in the A-AF phase. In the cubic lattice studied 
thus far, the C- and G-type orbital structures are degenerate, and it is unclear whether the 
orbital pattern in the x-y plane corresponds to the alternation of 3x 2 — r 2 and 3y 2 — r 2 
orbitals observed in experiments ||. To remedy the situation, some empirical facts observed 
in manganites become important: (i) The MnC>6 octahedra are slightly tilted from each other, 
leading to modifications in the hopping matrix. Among these modifications, the generation of 
a non-zero value for t* a is important, (ii) The lattice is not cubic, but the relation L m > L out 
holds, where L m = L x = L y and L out = L z . From experimental results |J, these numbers are 
estimated as L m = 4.12A and L ou t = 3.92A, indicating that the distortion with (^-symmetry 
occurs spontaneously. 

For the inclusion of the point (ii), in Eq. (|7|) is rewritten as 

Qsi = Qf + (1/V6)(25 ri - <5 xi - S yi ), (17) 

where Qg indicates the spontaneous distortion with (^-symmetry, given by 
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Figure 7: (a) Total energy as a function of i aa on the 2x2x2 lattice with L x = L y > L z for 
Jh = 8, A = 1.6 and J' = 0.05 (relaxation method). The solid and dashed curves denote the 
C- and G-orbital states, respectively, (b) Orbital structure in the A-AF phase for i aa = + 
and L x = L y > L z . These shapes are very close to purely 3x 2 — r 2 and 3y 2 — r 2 types. 

This length is rewritten in the non-dimensional form as = —r]X, where rj is a numerical 
factor given by rj = y/2/3{L- m — L out )/£jTj estimated as 0.5 by using the experimental data. 
Note that the hopping amplitude becomes different from those in the x-y plane due to this 
distortion. As discussed shortly in subsection 2.3, it is obtained as = (4t/3)(L in /L out ) 7 . As 
for J' along the z-direction, it is given by J'(L- m /L out ) 14: , since the superexchange interaction 
is proportional to the square of the hopping amplitude. 

Motivated by these observations, the energies for C- and G-type orbital structures were 
recalculated including this time a nonzero value for t| a in the magnetic A-AF state (see 
Fig. 7a)). In Fig. 8, the configuration of x 2 — y 2 orbitals is depicted. From this figure, it is 
intuitively understood that if the x 2 — y 2 orbitals in the upper and lower planes are tilted 
from the x-y plane as showing by arrows in the figure, there appears finite a hopping integral 
between adjacent x 2 — y 2 orbitals along the z-direction and the sign of this hopping amplitude 
is the same as that of i aa or t% a . Thus, in the real material, the tilting of the MnO"6 octahedra 
will always lead to a positive value for i aa and the results of Fig. 7(a) suggest that the C- 
type orbital structure should be stabilized in the real materials. The explicit shape of the 
occupied orbitals is shown in Fig. 7(b). The experimentally relevant C-type structure with 
the approximate alternation of 3x 2 — r 2 and 3y 2 — r 2 orbitals is indeed successfully obtained 
by this procedure. Although the octahedron tilting actually leads to a change of all hopping 
amplitudes, effect not including in this work, the present analysis is sufficient to show that the 
C-type orbital structure is stabilized in the A-AF magnetic phase when t| a is a small positive 
number, as it occurs in the real materials. Our investigations show that this mechanism to 
stabilize the C-type structure works also for the purely electronic model in the Hartree-Fock 
approximation . 
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Figure 8: Configuration of x 2 — y 2 orbitals connected by oxygen 2p orbitals in a 2x2x2 cluster. 
The open and hatched parts indicate the plus and minus sign in the orbitals, respectively. 
Arrows indicate the direction of the tilting. 

5 Discussion and Summary 

In this work, the Coulomb interaction term i? e i-ei has been neglected, but this detail needs 
further clarification. For this purpose, -ff e i- e i is written as 

flel-el = U ^i 7 T n i7l + U ' Y ™ia<J»ib<x' + J Y C ia CT C Lr' C ia<x' c ib CT) (19) 
i7 \aa' icra' 

where U is the intra-orbital Coulomb interaction, U' the inter-orbital Coulomb interaction, 
and J is the inter-orbital exchange interaction. For Mn-oxides, they are estimated as U = 
7eV, J = 2eV, and U' = 5eV [||, which are large compared to t. However, the result for the 
optimized distortion described in this paper, obtained without the Coulomb interactions, is 
not expected to change, since the energy gain due to the JT-distortion is maximized when 
a single e g -electron is present per site. This is essentially the same effect as produced by 
a short-range repulsion. In fact, the MC simulations show that the probability of double 
occupancy of a single orbital is negligible in the window of couplings where the A-type spin 
and C-type orbital state is stable. 

In order to confirm the above statement, the JT- and breathing-distortions were calculated 
as a function of U' by using the Exact Diagonalization method on a 2 x 2 cluster in the F-state 
in which U and J can be neglected. The result is shown in Fig. 9, where Qjt and Qbr are 
defined as 

QjT = (i/AO£\A& + Qii> (20) 

i 

and 

Q hl = (l/N)Y\Qu\, (21) 
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Figure 9: Mean-value of the breathing- and JT-mode distortions as a function of U' for the 
F-phase on a 2 x 2 cluster for Jh = 8 and A = 1.5. 



respectively. As expected, the mean value of the breathing-mode distortion is almost zero 
and only the JT-mode is active in the case of (3 = 2. It is noted that the dependence of Qjt 
on U' is very weak, indicating that the optimized distortion is not affected by the Coulomb 
interaction. The orbital arrangement in this 2x2 lattice is exactly the same as that in 
the x-y plane of the orbital structure for the A-AF phase in Fig. 5 and this arrangement is 
unchanged by the inclusion of U'. Note also that Qjt is gradually increased with the increase 
of U', although the dependence is weak. This suggests that the JT-distortion without U' is 
reproduced at smaller value of A if U' is included explicitly [16]. Thus, it is expected that 
the A-AF state will be stabilized at smaller A improving the comparison of our results with 
experiments. Based on all these observations, it is believed that the effect of the Coulomb 
interaction is not crucial for the appearance of the A-AF state with the proper orbital order. 
Another way to rationalize this result is that the integration of the JT-phonons at large A 
will likely induce Coulombic interactions dynamically. 

Finally, let us briefly discuss transitions induced by the application of external magnetic 
fields on undoped manganites. When the A-AF state is stabilized, the energy difference (per 
site) obtained in our study between the A-AF and F states is about i/100. As a consequence, 
magnetic fields of 20 ~ 50T could drive the transition from A-AF to F order accompanied by 
a change of orbital structure, interesting effect which may be observed in present magnetic 
field facilities. 

In summary, with the use of numerical techniques at n = 1, it has been shown that the A- 
AF state is stable in models with JT-phonons, using coupling values physically reasonable for 
LaMn03. Our results indicate that it is not necessary to include large Coulombic interactions 
in the calculations to reproduce experimental results for the manganites. Considering the 
small but important effect of the octahedra tilting of the real materials, the C-type orbital 
structure (with the alternation pattern of 2>x 2 — r 2 and 3y 2 — r 2 orbitals) has been successfully 
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reproduced for the A-AF phase in this context. 
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